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On the basis of the linear response transport theory, the general expressions for the thermoelec- 
tric transport coefficients, such as thermoelectric power (5), Nernst coefficient (y), and Thermal 
conductivity (k), are derived by using the Fermi liquid theory. The obtained expression is exact as 
for the most singular term in terms of 1/7^ (7k being the quasiparticle damping rate). We utilize 
the Ward identities for the heat velocity which is derived by the local energy conservation law. The 
derived expressions enable us to calculate various thermoelectric transport coefficients in a system- 
atic way, within the framework of the conserving approximation as Baym and Kadanoff. Thus, the 
present expressions are very useful for studying the strongly correlated electrons such as high-T c 
superconductors, organic metals, and heavy Fermion systems, where the current vertex correction 
(VC) is expected to play important roles. By using the derived expression, we calculate the thermal 
conductivity k in a free-dispersion model up to the second-order with respect to the on-site Coulomb 
potential U. We find that it is slightly enhanced due to the VC for the heat current, although the 
VC for electron current makes the conductivity (a) of this system diverge, reflecting the absence of 
the Umklapp process. 

PACS numbers: 74.25.Fy, 72.f0.-d, 72.10.Bg 



I. INTRODUCTION 

In general, the transport phenomena in metals are very important physical objects because they offer much in- 
formation on many-body electronic properties of the system. Especially in strongly correlated electron systems like 
high-Tc cuprates or heavy Fermion systems, various transport coefficients show striking non-Fermi liquid type be- 
haviors. Historically, theoretical studies on transport phenomena give considerable progresses in various fields of the 
condensed matter physics, such as Kondo problem and high-Tc superconductivity. 

According to the linear response theory [1-3] or the Kubo formula [4] , transport coefficients are given by correspond- 
ing current-current correlation functions. Thus, to study transport phenomena, we need to calculate the two-body 
green function with appropriate vertex corrections (VC's). Unfortunately, in many cases it is a difficult analytical or 
numerical work. Therefore, at the present stage, transport coefficients are usually studied within the relaxation time 
approximation (RTA), by dropping all the VC's. The effect of the VC can be included by the standard variational 
method by Ziman based on the Boltzmann transport theory [5] However, it is not so powerful for anisotropic corre- 
lated systems because there is no systematic way of choosing the trial function. Thus, it is desired to establish the 
microscopic transport theory based on the linear response formula. 

Based on the Kubo formula, Eliashberg derived a general expression for the dc-conductivity (a) in the Fermi liquid 
by taking the VC's into account, by performing the analytic continuation of the current-current correlation functions 
[6] . Based on the expression, Yamada and Yosida proved rigorously that a diverges even at finite temperatures if the 
Umklapp scattering process is absent. [7]. By generalizing the Eliashberg's theory including the outer magnetic field, 
the exact formulae for the Hall coefficient (Ra) [8] and the magnetoresistance (Ap/p) [9] in Fermi liquid systems were 
derived. By using these formulae, in principle, we can perform the conserving approximation for these coefficients 
with including appropriate VC's for currents [10]. 

In general, the VC's for currents are expected to be important especially in strongly correlated systems. For example, 
in high-Tc cuprates, the so-called Kohler's rule (|i?H| oc p° and Ap/p cx p~ 2 ) is strongly violated [11]. Moreover, 
i?H < in electron-doped compounds, although the shape of the Fermi surface is everywhere hole-like. These 
behaviors, which cannot be explained within the relaxation time approximation (RTA), had been an open problem in 
high-Tc cuprates. Based on the conserving approximation, we found that these anomalies are well reproduced by the 
VC's for electron currents [12-14]. The effect of the VC's, which are dropped in the RTA, becomes much important 
in a Fermi liquid with strong antiferromagnetic or superconducting fluctuations. 

However, as for the thermoelectric transport coefficients such as the thermoelectric power (TEP, S), the Nernst 
coefficient (v), and the thermal conductivity (k), we do not know useful expressions for analysis in the strongly 
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correlated systems so far. Here, the definition of S, k, and v under the magnetic filed B parallel to the z-axis are 
given by 

S = —E x j d x T, 

« = -Q x /d x T, (1) 
v = -E y /Bd x T, 

where Q is the heat current. Unfortunately, the conserving approximation for these coefficients is not practicable 
because we do not know how to calculate the VC's for them. Thus, at the present stage, only the RTA is accessible, 
which will be insufficient for a reliable analysis of the strongly correlated systems because VC's should be included. 

There is a long history of the microscopic study on thermoelectric transport phenomena. To this problem, we 
cannot apply the Kubo formula to the electronic conductivity naively because there is no Hamiltonian which describes 
a temperature gradient d x T . In 1964 Luttinger gave a microscopic proof that thermoelectric transport coefficients 
are given by the corresponding current-current correlation function [1]. Later, Mahan et al. much simplified the 
Luttinger's expression in the case of electron-phonon and electron-impurity interactions [3]. However, the analysis 
of the VC for the heat current for electron-electron interacting systems is still an open problem, which is necessary 
to go beyond the RTA. This analysis will be more complicated and profound than that for the electric conductivity 
performed by Eliashberg [6] . 

In the present paper, we derive the thermoelectric transport coefficients by performing the analytic continuations of 
the current-current correlation functions, on the basis of the linear response theory developed by Luttinger or Mahan. 
Our expressions are valid for general two-body interactions. The VC for the heat current is given without ambiguity 
by the Ward identity with respect to the local energy conservation law. The derived expressions are "exact" as for 
the most divergent term with respect to 7j7 , where 7k is the quasiparticle damping rate. The present work enables 
us to perform the "conserving approximation" for S, v and k, which is highly demanded to avoid unphysical results. 
Actually, the VC's would totally modify the behavior of these quantities in strongly correlated electron systems, as it 
does for the Hall effect and the magnetoresistance. In this respect, the RTA is unsatisfactory because all the current 
VC's are neglected there. 

We note that Langer studied the Ward identity for the heat current, and discussed the thermal conductivity [15]. 
However, the derived Ward identity was not correct as explained in §111, although it did not influence the thermal 
conductivity at lower temperatures fortunately. We derive the correct Ward identity in §111, and give expressions for 
k, S and v. 

In thermodynamics, the TEP of metals becomes zero at absolute zero temperature, which is the consequence of 
"the third law of the thermodynamics" . (As is well known, the third law also tells that the heat capacity vanishes at 
T = 0.) Similarly, both v and k also become zero at T = if the quasiparticle relaxation time r (= I/27) is finite 
at T = due to impurity scatterings. Unfortunately, these indisputable facts are nontrivial in a naive perturbation 
study once the electron-electron correlations are set in. In the present work, we derive the general expression for S 
which automatically satisfy S(T = 0) = owing to the Ward identity for the heat velocity. 

In high-Tc cuprates, the Nernst coefficient v increases drastically below the pseudo-gap temperature, which is never 
possible to explain within the RTA [16]. According to recent theoretical works, superconducting (SC) fluctuation is 
one of the promising origins of the pseudo-gap phenomena [17-20]. Based on the opinion, we studied v for high-Tc 
cuprates using the general expression derived in the present paper [21]. Then, we could reproduce the rapid increase 
of v only when the VC's due to the strong antifcrromagnetic and superconducting fluctuations are taken into account. 
This work strongly suggests that the origin of the pseudo-gap phenomena in high-Tc cuprates is the strong d-wave 
superconducting fluctuations. 

In the case of heavy Fermion systems, the TEP takes an enhanced value around the coherent temperature, and its 
sign changes in some compounds at lower temperatures. Such an interesting non-Fermi liquid like behavior is mainly 
attributed to a huge energy dependence of the relaxation time, r(e) (= l/27(e)), due to the Kondo resonance. This 
phenomenon was studied by using the dynamical mean field theory [22,23]. Also, the TEP in the Kondo insulator 
was studied in rcf. [24] in detail. 

The contents of this paper are as follows: In §11, we develop the linear response theory for thermoelectric transport 
coefficients. By performing the analytic continuation, we derive the general formula of S and k in the presence of 
the on-site Coulomb potential U. In §111, we derive the Ward identity for the heat velocity which is valid for general 
two-body interactions, by using the local energy conservation law. The Ward identity assures that the expressions for 
S and k derived in the previous section are valid even if the range of the interaction is finite, as for the most divergent 
terms with respect to 7 _1 . In §IV, the general formula for v is derived. It is rather a complicated task because the 
gauge invariance should be maintained. In §V, we calculate k in a spherical correlated electron system in the absence 
of Umklapp process, and obtain its exact expression by including the VC's within the second order perturbation. The 
physical meaning of the VC is discussed. Finally, the summary of the present work is shortly expressed in §VI. 
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II. LINEAR RESPONSE THEORY FOR S AND k 



First, we shortly summarize the linear response theory for thermoelectric transport coefficients, initiated by Lut- 
tingcr. Here we consider the situation that both the electron current J 1 = J and the heat current S 2 = Q are caused 
by the external forces X 1 = E/T and X 2 = V (7^), where E is the electric field. In the linear response, 

ji = J2 L lm (B)X m , (2) 

m=l,2 

where l,m= 1,2. Because the relation ^f- = J 1 ■ X 1 (S e being the entropy) is satisfied in the present definition, 
the tensor L lm (B) satisfy the Onsager relation; L l ™(Ti) = L™^(— B), where /j,, v = x, y, z [25]. 

According to the quantum mechanics, the electron current operator j and the heat current one 2 are given by [2] 

j(r,) =i[H,ep(r l )r t ], (3) 
j Q (r i )=j E (r i )-^j(r i ), (4) 
] E {v i )=i[H,h{v i )v l l (5) 



where e(< 0) is the charge of an electron, H is the Hamiltonian without external fields X 1 , p(r^) = cJ.(rj)c (T (rj) is 
the density operator (a being the spin suffix), and h(ri) is the local Hamiltonian by which H is given by H — h(ri). 

By using these current operators, J and Q are given by J(ri) = O'(rO) an d Q{ r i) — (j®( r i)), respectively. 

To derive the expressions for various conductivities microscopically, we introduce the virtual external potential term 
F which causes the currents X 1 (I = 1,2). Then, the "total Hamiltonian" is expressed as Ht = H + F ■ e (-* w + l5 ) t j 
where 6 > is a infinitesimally small constant. According to the linear response theory [1,4], the current J 1 at t is 
given by 



= -if dt'(g i (c l = 0,t),F(t')})e(- i " (6) 

dF OS — — 

Because of the relation — — = T—^- = T^^^f ' * ne expression for L lm is given by [2] 



-Iva rim 



L tm (cu + id) , (7) 

U! = 



-T 







l!™{k*t) = — / (T T ^(q = 0, r)j™(q = 0, r = 0)) , (8) 

^1 Jo 

where (3 — 1/T and [i, v = x,y. T T is the r-ordering operator, and u>i — 2irTl (I being the integer) is the bosonic 
Matsubara frequency. By writing the diagonal component of L lm as L lm , a, S and k are given by [2] 

- = fU" (9) 
S=i^T = ^ — , (10) 



eTL ll T 2 (1 

1 

2^2 



1/ 22 _IW, 



where e (e < 0) is the charge of an electron. 

Hereafter, we analyze the function L lm {iuji) given by eq.(8) at first, and perform the analytic continuation to derive 
L lm by eq.(7). We study a tight-binding model with two-body interactions, which is expressed in the absence of the 
magnetic field as: 
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H = H Q + H int , (12) 
tfo=E £ kcL^, (13) 



k,<r 

= 2 H l7 <7<T'( < l) c k+q,<r C k'-q,-<7^k' -<rCk,<7- ( 14 ) 
kk'qerer' 

In eq.(13), e k = ii.o e * k '( ri_r °), where ijj is the hopping parameter between and Yj. U aa i(q) represents the 
electron-electron correlation between a and a' spins. For example, U aa >(q) = U8 a - a ' for the on-site Coulomb 
interaction. 

The one- particle Green function is given by Gk(e) = l/(e + [i — e k — Sk(e)), where Sk(e) is the self-energy and 
ji is the chemical potential. In a Fermi liquid, 7k <C T is satisfied at sufficiently low temperatures because of the 
relation 7k oc T 2 [26]. In such a temperature region, the following quasiparticle representation of the Green function 
is possible: 

GkM = ^ -, (15) 

where = (1 — ^jEk(e)) -1 is the renormalization factor, e k = Zk(ek + ^k(O) — /u), 7 k = Zk7k and 7k = Im£k(— «0), 
respectively. 

According to eq.(3), the electron current operator for eq.(12) is given by 

J(p) = e ^ ^c k -p/2, ff c k+ P /2, ff , (16) 

k,cr 

where = Vfee k . Apparently, j(p) is a one-body operator. 

In the same way, we consider the heat current operator defined by cq.(12): In the case of the on-site Coulomb 
interaction, for simplicity, it is obtained after a long but straightforward calculation as 



k,<T 

+ ^ J2 + W t-q/2) C L q /2,/k+q/2, ff i + q /2,-/t-q/2,- ff , (17) 

kk'qer 

which contains a two-body term in the case of U ^ 0. It becomes more complicated for general long-range potential. 
This fact seems to make the analysis of the thermoelectric coefficient very difficult. 

Fortunately, as shown in Appendix A, eq.(17) can be transformed into the following simple one-body operator form 
by using the kinetic equation: 

J Q (P = 0,«l) = £/^ r }™l - -jp) 4>K ct (t') 

= T *( e ™ + wi/2)i$ • 4 ;CT (e n )c k;CT (e„ + w,), (18) 

where and e„ are boson and fermion Matsubara frequencies, respectively. The case of the nonlocal electron-electron 
interaction is discussed in the next section by constructing the Ward identity for the heat velocity. 
By using eq.(18), we obtain the expression for L 12 without the magnetic field as follows: 



T 2 e v-^ 1 

L ll( iuj i) = 2_/( e ™ + o w O w k ([i fl , k(en,wOA kl/ (e„,a; i ), (19) 

k,e„ 

A kl ,(e n ,w;) = v^ v +T ^ r(ke„,k'e„/;u;()gk'(en',^) w kv> (20) 



k',e„ 



where gk(e n ,uji) = Gk(«£ n -\- iu)i)G\;_{ie n ). Ak„(e„,o;;) and r(ke„, k'e„/ ; u>i) are the three- and four-point vertices 
respectively, which are expressed in Fig. 1. They are reducible with respect to the particle-hole channel. Note that 
we put the outer momentum p = in cq.(20) because we are interested in the dc-conductivity. 



4 



k,e k,e+co 

= r(ke,k'e';co) 



k'£' k'£'+CO 
k,E k,E+CO 

= A kv (e;w) 



FIG. 1. The four-point vertex correction T and the three-point vertex function A„, respectively. 



The expression (19) for L 12 derived for the on-site Coulomb interaction is equal to that for a system with the 
impurity scattering and the electron-phonon interaction derived by Jonson and Mahan. [3]. In the next section, we 
will show that the expression is also valid for general types of two-body interactions as for the most divergent term 
with respect to 7 _1 on the basis of the Ward identity. 

The dc-TEP is obtained by the analytic continuation of eq.(19) with respect to iu>i, by taking all the VC's into 
account. (The analysis on the VC's in ref. [3] is insufficient.) In the present work, we perform the analytic continuation 
rigorously by referring to the Eliashberg's procedure in ref. [6]. Next, by using the Ward identity (26), we derive the 
simple expression for the TEP within the most divergent term with respect to 7 . 

After the analytic continuation of eq.(19), L 12 v (+i8) of order (9(7 _1 ) is given by 

Lll(+iS) = eTj2 J * (-§Q q^(e)\G k (e)\ 2 JUe), (21) 
where the total electron current with VC's and the quasiparticle heat velocity q^ v are respectively given by 

Me) = Mt) + J2 J ^(-^)r 22 (k e ,kV)^( e 'H'(0, (22) 

Me)=€+ E /^T*(ke,kV) 5 £VK'> (23) 
<fk(e)=«£+ £ /^T 2l (k £ ,kV) 5 «( e ')^, (24) 

k'.i=l,3 

where g^\ e ) — {^(e)} 2 , g^\e) = |Gk(e)| 2 and 5 k ^( e ) = {'-■'kH 6 )} 2 ! respectively. The definition of the four-point 
vertex T m,z (pe, p'e') is given in ref. [6], which are listed in Appendix B. In general, T 2t is well approximated at lower 
temperatures as (T u + T 3 *)/2 [6]. Thus, taking the Ward identity for electron current is taken into account [26,27], 
Wk(e) is simply given by 

#k(e) = V fe (e k + ReE k (e)). (25) 

Next, we consider <?k(e) defined in eq.(24): By seeing its functional form, the relation qy(e = 0) =0 is nontrivial. 
However, if q\(e = 0) were nonzero, then L 12 in eq.(19) would be proportional to T-y^ 1 . In this case, S = eL 12 /crT 2 
diverges at T = 0, which contradicts "the third law of the thermodynamics"! In this sense, eq.(24) is too primitive 
for a reliable (numerical) analysis at lower temperatures. 

Fortunately, by noticing that T 2% = (T l1 +T 3l )/2 at lower temperatures, the quasiparticle heat velocity <fk(e) given 
in eq.(24) can be rewritten in a simple form as 

?k(e) = «T k (e), (26) 

where Vk(e) is given in eq.(25). Equation (26) is the Ward identity which will be derived from the local energy 
conservation law in the next section. This Ward identity leads to L 12 oc T 3 7~ 1 because of q\(e = 0) = 0, so the 
difficulty in analyzing the TEP towards T — > is removed. 

In the same way, we derive the exact formula for the thermal conductivity within the most divergent term with 
respect to 7 _1 . By the similar way to the derivation of eq.(19), we obtain that 
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T 1 

LfA^i) = —J2 i{ ~ €n + oWi)<,flk(cn,wi)A£,(e n ,wi), (27) 
A kl( e ™> w O = *( e « + \^ML + T X! r ( ke n,k'e„,;a);).9k'(en',^) • i{e n > + \^i)v^, v . (28) 



k' 

After the analytic continuation, we find that 



Lf v (+iS) =Tj2j ^(~%) ^( e )|G k ( £ )| 2 Q k „(6), (29) 
Qk(e)=? k (e) + X: /^(-g)r 22 (k e ,kV)^( e ')^( £ '), (30) 

k' ^ ' 

where Qfc(e) is the total heat current with VC's, which is given by the analytic continuation of A^ (ie„; iu>i) from the 

region e„ + uii > and e„ < 0, and by taking the limit w — > 0. We stress that Qk(e = 0) = at zero temperature, 
because the e'-integration range in cq.(30) is restricted to within |e'| <~ T due to the thermal factors in T 22 given by 
the analytic continuation; see Appendix B. This fact leads to the result that L 22 v (+iS) given by eq.(29) is proportional 
to T 3 7~ 1 as T — ► 0. As a result, at lower temperatures, the relation /t oc T-f^ 1 is assured in the present analysis. 
In summary, (i) the TEP is given by S — —r—L 12 , where L 12 is given in eq.(21) and J, v and q are given by eq.(22), 

L 22 

eq.(25) and eq.(26), respectively, (ii) the thermal conductivity is given by n = — ^ — TS 2 a, where L 22 is given in 

eq.(29). At lower temperatures, the first term L 22 /T 2 is dominant because it is proportional to Tj^ 1 whereas the 
second term is proportional to T 3 7 _1 . 

Finally, we note that the conductivity <r is given by [6] 

<V = e 2 Y,j f (-%) ^( £ )|Gk(£)| 2 Jk,( £ ), (31) 

k 

where Jk(e) is given by eq.(22). 



III. GENERALIZED WARD IDENTITY 



As we discussed in §1, the TEP becomes zero at zero temperature, which are ensured by "the third law of the 
thermodynamics" . This fact means that L 12 given by eq.(21) should be proportional to T 3 7 _1 as T — > 0. This relation 
is ensured if the Ward identity (26) is satisfied exactly. In this section, we derive the generalized Ward identity for 
general types of the two-body interactions, by noticing the local energy conservation law; J^/i + V • = 0, where h(z) 
is the local part of the Hamiltonian. The obtained heat velocity has a correction term Aq (see eq.(38)), which turn 
out to be negligible for transport coefficients at lower temperatures. The present derivation is analogous to the proof 
for the generalized Ward identity for the electron current which is described in ref. [28]. In this section, we drop the 
spin suffixes for simplicity of the description. In the same reason, we put \i = because it will not cause a confusion. 

Here, we introduce the four-dimensional heat velocity; (h(z),j®(z)) = (jo(z), ji(z), J?( z ))- Then, we 

consider the following function X® (^ = 0^3): 

X${x,y,z) = (T T j$(z)c{x)c\y)) 

= JI G(x, x')AQ(x', y>, z)G(y>, y)dVdV , (32) 

where Ajj (a/, y' , z) is the three-point vertex function with respect to the heat velocity. For the simplicity of the descrip- 
tion, we assume hereafter that x, y and z are continuous variables, not discrete ones. Because of the translationally 
invariance of the system, we can write as 

X$(x, y,z)= f f X®(k+ p, k)e l(k(x -y*> + P( x - z »d 4 kd 4 p. (33) 



G 



k,e k,e 



k,e k,e 



k,e k,e 



+ 



1 Talk's') 



%(£) 



V 



s'v k °. 



>ev k u 

= 8Vk(e k ^£ k (8)) 

FIG. 2. The Ward identity for the heat velocity <fk(e) derived in this section. This identity assures that g k (e = 0) = 0. 



Here, we use the following four-dimensional notations for x, p and k; (t,xi,x 2 ,x 3 ) = x^, (— cj,pi,p 2 ,P3) = Pp, and 
(-e,ki,k 2 ,h) = kf,. 

For the moment, we assume that h(z) is a local operator i.e., the two-body interaction f/(x — y) is a (5-function 
type. (This restriction on h(z) is released later in the present section.) Because the relation i[h(z), c(x)]S(zq — xq) — 
g|^c(x)<5 4 (z — x) is satisfied, then 



v=0 



dz v 



dz u 



c(x)c t (y) 



+ (T r [/i(z),c(x)]c t (y))<5(z - x ) 
+ (T r c(x)[/i(z), C t(y)]) ( 5(z -y ) 

- y) • <5 4 (z - z) + A G(a; _ y ) . 5 \y - z), 
oxq oyo 



(34) 



where <5 4 (a;) = S(x)S(xo). In the transformation, the local energy conservation law is taken into account. Performing 
the Fourier transformation of eq.(34), we get 



5>„A«(fc +p, k) = (el + £(fc))(fco + Po) - + +p))Ao. 



(35) 



i/=0 



By putting p M = for /x 7^ i (i =1,2 or 3) in eq.(35) and taking the limit pi — > 0, we obtain the z-component of the 
heat velocity qi(k) as follows: 



gi(fc) = lim Af(k+pi,k) 

Pi— >0 



(36) 



for z = 1,2,3, which is equivalent to the Ward identity for the heat velocity, eq.(26). By constructing the Bethe- 
Salpeter equation [26,27], q(k, e) in eq.(36) is expressed by using the fc-limit four-point vertex r fe (ke,k'e') as follows: 



&(e)=cV fc (c£ + S(p)) 

= ^ + pE/ ^^ ( ke,k'e'){G k , ( e' ) 2 } fe • e'tfc 



(37) 



in terms of the zero-temperature perturbation theory, which is diagrammatically shown in fig. 2. In the finite temper- 
ature perturbation theory, q p (e) is expressed as eq.(24). We stress that eq.(37) is satisfied only when we take account 
of all the diagrams for T k (ke, k'e') which are given by the functional derivative ST,/5G. (see Appendix D.) 
Finally, we discuss the following two restrictions assumed in the proof of the Ward identity: 

(i) In the above discussion, we assumed that the potential term in h(z) is local, which is not true if the range of the 
interaction is finite. In this case, a correction term C(x,y;z) given in eq.(C9) should be added to (34), as discussed 
in Appendix C. This correction term gives rise to the additional heat velocity Aq(k) given by eq.(C13). Note that 
A<fk(e) vanishes identically in the case of the on-site Coulomb potential. As a result, the Ward identity for the heat 
velocity in the case of general two-body interactions is given by 
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q(k, e) = ev k (e) + A<f k (e). 



(38) 



Fortunately, Aq\k) does not contribute to the transport coefficients as discussed in Appendix C. As a result, we 
can use the expression for q(k) in eq.(37) for the purpose of calculating S, k and u, even if the potential U(k) is 
momentum-dependent . 

(ii) Here we treated the space variables like x, y and z as continuous ones for the simplicity. However, it is easy to 
perform the similar analysis for the tight-binding model, by replacing the derivative of x with the differentiation. For 
example, the local energy conservation law is expressed as §ih(ri) + (j(r; + i) — j(r/_i))/2a = 0, where a is the lattice 
spacing. We stress that the Ward identity in eq.(36) is rigorous also in the case of the tight-binding model. Note that 
we give the another proof for the Ward identity based on the diagrammatic technique in Appendix D, which is valid 
for general tight-binding models. 

We comment that Langer studied the Ward identity for the heat velocity in ref. [15]. Unfortunately, because of a 
mistake, an extra factor — i7kG k _1 (w) should be subtracted from the r.h.s. of his Ward identity (i.e., the heat velocity), 
eq.(3.31) in ref. [15]; see eq.(37) or (38) in the present paper. In fact, the factor k • (k + q) in cqs. (3.16) and (3.22) 
(in eqs. (3.19) and (3.23)) of ref. [15] should be replaced with k 2 (with (k + q) 2 ). This failure, which is fortunately 
not serious in studying the transport coefficients at lower temperatures, becomes manifest if one studies the Ward 
identity in terms of the x-representation, like in the present study. 



IV. FORMULA FOR THE NERNST COEFFICIENT 



According to the linear response theory, the Nernst coefficient is given by 

/B, (39) 



^-Stanfe 



where a(B) = L 12 (B)/T 2 is called the Peltier tensor (J = a(— VT)), and tan#n = &xy/<J is the Hall angle. Note 
that a yx (B) = L y l(B)/T 2 = L 2 ^(-B)/T 2 = -L 21 y (B)/T 2 because of the Onsager relation. In addition, a yx (B) = 
— a xy (B) in the presence of the four- fold symmetry along the magnetic filed B. 

In this section, we investigate the off-diagonal Peltier coefficient a yx due to the Lorentz force to derive the expression 
for the Nernst coefficient. Up to now, the general expression for the Hall coefficient [8,29] and that for the magne- 
toresistance [9] were derived by using the Fermi liquid theory based on the Kubo formula. These works enabled us to 
perform numerical calculations for the Hubbard model [12,14,31] within the conservation approximation as Baym and 
Kadanoff [10]. Hereafter, we derive the general expression for the Nernst coefficient by using the technique developed 
in refs. [8], [9] and [29]. 

For the present purpose, we have to include the external magnetic field. In the presence of the vector potential, the 
hopping parameter tij in eq.(13) is multiplied by the Peierls phase factor [9]: 

Uj -> Uj exp[ie(Ai + Aj) ■ (r* - rj)/2], (40) 

where Aj is the external vector potential at r i} and e(< 0) is the charge of an electron, Here we introduce A, as 

A 4 = Ae 4 P- r % (41) 

where A is a constant vector. In this case, the magnetic field is given by B = ip x A in the uniform limit, i.e., |p| -C 1 
[8,9,29]. Bearing eqs. (40) and (41) in mind, the current operator defined by eq.(3) and the Hamiltonian are given by 
[9] 

3„(P = 0) = J AO) - eA a ■ J„a(-P), (42) 
H B = H B=0 - eA a j a (-p), (43) 

in the tight-binding model up to O(A). Here and hereafter, the summation with respect to the suffix which appears 
twice is taken implicitly. j a (p) is given in eq.(16), and 

j a p(p) = e^( 9 « 9 /3 £ k)c k -p/2 c k+p/2- (44) 

k 

To derive the Nernst coefficient, we have to calculate the L\} under the magnetic field, which is given by 
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FIG. 3. 

All the diagrams for a pv . The symbol "op" on each line represents the momentum derivative, "d p " . The notations 

in the diagrams are explained in ref. [9] in detail. 



L^(p,^;A) = -Z^ /3 dre-^iTr{e-^T r j« B (p,r)^(0,0)}, (45) 

and take derivative of eq.(45) with respect to q p and A a up to the first order. Z is the partition function. By taking 
eqs.(42) and (43) into account, the A a , ^-derivative of L^(p,iw/; A) is given by 

-a^iw,) = J2 C%(iLJi)(i Pp ■ A a ), (46) 
op 

-i d 2 



T^d^dA a L ^ P ^ A) 



(47) 

p=A=0 



Below, we see that the dc-Peltier coefficient a yx is given by the analytic continuation of a yx (iu>i). 

The diagrammatic expression for Cgy(iu)i) is very complicated, containing six-point vertices. Fortunately, as for 
the most divergent term with respect to 7 _1 , they can be collected into a small number of simpler diagrams as shown 
in fig. 3 by taking the Ward identity into account [8,9,29]. We can perform the present calculation for Cgy(iwi) in 
a similar way to that for o xy in ref. [8] , only by replacing j p ,= x 

with jQ =x and using the Ward idntity for the heat 

velocity. As a result, we obtain the following result: 

C%{iui) = ^E A k*(w^) {[G o p G + ]d a [G da G+ \ d p] Ak^e„^) 

k,ie n 

+— E ^ x (ien,iui)[d p G + ■ d a G - d a G+ ■ d p G]k± v {ie n ;iuJi), (48) 

where G + = Gk(e; + uji), G = Gk(e;), and [A Q a B] = A-d a B — B-d a A. A® x (ie n ; iu>i) is the three point vertex for 
the heat current given in eq.(28). Equation (48) is described in (a)-(d) fig. 3. Here, we neglect the diagrams (e) and 
(f) because their contribution is less singular with respect to 7 _1 [8]. Here, we assume that the magnetic field B is 
parallel to the z-axis. Then, we can easily check that C^ a y (i<^i) in eq.(48) is expressed as CP y =x ' ,J=v (iivi) ■ e zprT . This 
fact assures that a xy (iuJi) given in eq.(47) is gauge-invariant, that is, 

—a yx (iuji) = C^y X,cr=y (i(Ji) ■ B z . (49) 
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We note that eq.(48) is equivalent to a xy {iu>i) /T if A® y is replaced with Akj,; see "A" and "B" in p. 632 of ref. [8]. 

In performing the analytic continuation of eq.(48), the most divergent term with resect to 7 -1 is given by the 
replacements G(e n +uji) — > G R and G(e„) — ► G A . Taking account of the relation d p G R (e) — G R (e) 2 (t>kp(e) — «<9 P 7k(e)), 
the dc-Peltier coefficient a xy (uj + i6)\u=o is obtained as 



B - ^E/ % (-%) [\^G k (e)\\G k (e)\ 2 Q kx [v kx d y - v ky d x ]J ky 
+ |G k (e)| 4 Q kx Jkj / [-Wk 2; 9 y 7k + UkjA7k] j, (50) 

where Q k (e) is the total heat current introduced in eq.(30) We stress that Q k (e = 0) = at T = 0, as is discussed in 
§11. 

It is instructive to make a comparison between a xy and <r xy /eT: The latter is given by eq.(50) by replacing Q x 
with J x . In this case, the second term of eq.(50), which contains the /c-dcrivativc of 7k, vanishes identically because 
of the Onsager relation a xy (H) = — a yx (B). As a result, the general expression for a xy given in eq.(3.38) of ref. [8] is 
reproduced. 

If the system has the four-fold symmetry along the z-axis, then a xy (B) = — a yx (B). In this case, considering that 
|ImG k (e)||G k (e)| 2 = |G k (e)| 4 7k, eq.(50) can be rewritten as [12,13] 

-a yx = a xv = B.'^Y i J |L |ImGk( e )||G k ( e )| 2 7k ( e )A k ( e ), (51) 

A k (e) = (Q k (e) x (v k (e) x V) z ( J k (e)/ 7k (e)) ) 

= \Me)\± (&(e) x A (J k (e)/ 7 k(e))) , (52) 



where w k |± = \/v kx + v^_ x , and fc|| is the momentum on the xy-planc along the Fermi surface, i.e., along the vector 
ej| = (e z x Vk)/\vk\±- As noted above, eq.(51) becomes <r xy /eT by replacing Q x with J x ; see eq.(22) in ref. [12]. 
It is notable that A k (e) in eq.(52) is rewritten as 

qqJ , _ _ n g 

(7k(e)/|«U)^k(e) = (QkxJk. + QkyJky)-g± + [Qu x J k J z ^log(|J k |/7k), (53) 

where 8^. = tan _1 (Jk2,/Jky). In an interacting system without rotational symmetry, the second term with ku -derivative 
of 7k does not vanish in general since J is not parallel to Q owing to the VC's by T 22 . In contrast, Wk(e) = <fk(e)/e 
because of the Ward identity. In ref. [21], based on the fluctuation-exchange (FLEX) +T- matrix approximation, we 
studied the Nernst coefficient of the square lattice Hubbard model as an effective model for high-T c cuprates. Wc 
found that the second term of eq.(53) gives the huge contribution in the pseudo-gap region if the VC's for currents 
are taken into account in a conserving way. As a result, the origin of the abrupt increase of the Nernst coefficient 
under the pseudo-gap temperature is well understood. 



V. DISCUSSIONS 



A. Vertex Correction for Thermal Conductivity 

In previous sections, we studied various analytical properties for qu.(e) or Q k (e), using the Ward identity for the 
heat velocity derived in §111. In this subsection, wc study a free dispersion model (ek = k 2 /2m) in the presence of the 
electron-electron interaction without Umklapp processes. This situation will be realized in a tight-binding Hubbard 
model when the density of carrier is low; n«l. Here, we explicitly calculate the total heat current Qk(e) in terms 
of the conserving approximation. The present result explicitly shows that Jk(e) 7^ Q k {e)/e. 

Next, as a useful application of the expression for the transport coefficients derived in previous sections, we study 
the thermal conductivity k in a free-dispersion model. Because of the absence of Umklapp processes, the (T 2 -term 
of the) resistivity p of this system should be zero even at finite temperatures. In a microscopic study based on the 
Kubo formula, this physical requirement is recovered by taking account of all the VC's for the current given by the 
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(a) 



(b) 



(c) 




FIG. 4. The vertex corrections (by 722) for the heat/electron current in the second order perturbation theory. 



Ward identity [7] . On the other hand, the thermal conductivity is finite even in the absence of the Umklapp processes 
because heat currents are not conserved in the elastic normal scattering processes. Hereafter, we derive the Tj^ 1 - 
linear term of n in the free dispersion model in terms of the conserving approximation. For this purpose, we can 
drop the second term of eq.(ll) because L 12 L 21 /T 2 L n = TS 2 a ~ 0(T 3 7 1 ). The obtained result is exact within the 
second order perturbation with respect to U. 

First we consider the second order VC's as shown in fig. 4. Because Qk(e = 0) = 0, we can write Qk(e) = <?k • e up 

->(a-c) 

to O(e). The correction terms given by (a-c) in fig. 4, AQ k (e), are given by 



AQ k r) ( e )^ 2 £/ 



4 



, e ' - e , e' 
cth— — - th— 
2T 2T 



Ti%(e,e')\G k ,(e')\ 2 Q k/ (e') 



(54) 



where r — a,b, c. T^,(e, e') is a VC which is classified as T 22 : Their functional form are given by 



Ti%(e,e') = -Im X i\,(e-e') 



, u + e , w + e' 
-th— — +th- 



2T 



a; + e 



2T 



Pk+p(e + ^)Pk'+ P (e' + w), 



2T 



2T 



Pk+ P (e + w)p k '+p(e' - w), 



(55) 

(56) 
(57) 



where pk(e) = -ImG k (e — i<5) and Gk(e)| 2 = 7rpk(e)/7k(e)- By expanding eq.(54) with respect to e and T up to 



9 /" e' — e e' 
0(e 2 ,T 2 ) as was discussed in rcf. [7] and noticing that — / de'[cth— — ^^K 6 ' ~ e -' e ' 

that 



-{ttT) 2 , we obtain 



e=0 



AQ k a) = TT^ 2 E^+p(0)^ +P (0)^(0) 



k' P 

A<& 6) = A<&°\ 



e 2 + (ttT) s 
27k' (e) 



Ck', 



(58) 

(59) 
(60) 



In deriving eq.(60), we have changed the integration variables (k',e') — > (— k', — e'), and used the relation p- k (0) = 
p k (0) and Q_k(— e) = Qk(e)- In general, within the FLEX approximation, the Aslamazov-Larkin (AL) type VC's by 
T 22 , which correspond to (b) and (c), turn out to cancel out for the heat current. 
In the same way, the imaginary part of the self-energy, 7k(e), is given by 



U 



7k(e) = (e 2 + (ttT) 2 )— ^7rp k+p (0)p k , + p(0)p k ,(0), 



(61) 



k'p 



which is shown in fig. 5. 

In a spherical system, we can put Ck = C jj^ on the Fermi surface. Then, the total correction for Q is given by 
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E k (e) = !^ 
: — > 



FIG. 5. The self-energy given by the second order perturbation. 



AQ k = ]T AQ 



(r) 
k 



= 3^E-°k+p(0W +P (0W(0) -C^-, (62) 

k'p 

Z = ^Pk + p(0W +P (0W(0). (63) 
k'p 

Here, we put k = (0, 0, k F ). Then, the z-component of eq.(62) is given by 

j dk'd6 kl d(f> kl J dpd9 p d(f>pk' 2 p 2 sin9 k > sin9 p ■ 5(ev - ^)5(e k+p - ^)5(e k/+p - ^) ■ cos8 k > . (64) 

Note that in a free dispersion model, p k (0) = z5{z{e\ 5 _ — fj,)) = 6(6^ — fi), where e k — k 2 /2m and z is the renormalization 
factor. By performing fc'-integration, ^-integration, and <f) p , <f) k > -integrations successively, eq.(64) becomes 



4TTZ 3 m 3 e f in f 2k r sin ^/(i-cose,,) sin fl fc / cos ^ 

— ^ — / dBu / dp 



Ak 2 F sm 2 By -2p 2 (l-cos6> fe /) 



= Anrrfe J* sin ^ cos ^ = 8^ 



In the same way, Z is calculated as 



7*' 

Jo 



As a result, AQk is given by 



o >/2(l - cos6» fe /) 
8ttto 3 . (66) 



By solving the Bethe-Salpeter equation, Q = <f + AQ, we get 

Qk = (68) 

where <fk = ev k . As a result, the thermal conductivity within the second-order perturbation theory is given by 

k = (69) 

k° = , 70 

where k° is the result of the RTA, where VC's are neglected. Note that m in eq.(70) is unrenormalized, and n is the 
number of electrons in a unit volume. Finally, performing the momentum summations in eq.(61), 7 of order U 2 is 
given by 
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(71) 



In conclusion, the vertex corrections slightly enhances (by | times) the thermal conductivity in a three-dimensional 
frcc-dispcrsion model within the second order perturbation theory. 

It is instructive to make a comparison between the role of VC's for the heat current and that for the electron 
current. The VC's for the electron current which correspond to fig. 4 (a)-(c) are given by 

A j£° = A j£ 6) = Aj£ c) 

= t/ 2 ]T 7r Pk+p (0W +p (0W(0) f2 ± { f/ J v , (72) 

which was already derived in ref. [7]. Here we put Jk = D-^ on the Fermi surface. Performing all the momentum 
integrations in the spherical case as before, we find that 



A J k = AJ ) 



a,b,c 



= | E Pk+ P (0W + p(0W (o) • D^- 

k'p 

- D b (73) 

where Z is given in eq.(63). As a result, the solution of the Bethe-Salpeter equation J = v + A J is given by J = oo, 
which means that the conductivity diverges in the absence of the Umklapp processes, even at finite temperatures. 
Thus, the important result in ref. [7] is recovered. On the other hand, the thermal conductivity does not diverge even 
in the absence of the Umklapp processes, because the normal scattering process attenuates the heat current. 



B. The TEP and the Nernst Coefficient 



In this section, we discuss the effect of the anisotropy as well as the role of the VC's for the TEP and the Nernst 
coefficient. First, we discuss the validity of the Mott formula for S [30] which is given by 



3e 



d\na(E) 
dE 



(74) 



It is easy to see that eq.(74) is valid even in the presence of Coulomb interactions, if we define a(e) = 
&2 J2k |Ck(e)| 2 Vka;(e) Jka;(e) [3]: a and S given by eqs.(21) and (31) are rewritten using a(e) as 



-/?(- 



de 



S = 



7T 



de 



ecr(e), 



At sufficiently lower temperatures, eqs.(75) and (76) become 

a = a(0), 

_ _ 7r 2 fc|T da(e) 



3ea 



de 



(75) 
(76) 

(77) 
(78) 



e=0 



As a result, Mott formula is also satisfied in the case of electron-electron interaction. Note that the renormalization 
factor z does not appear in cq.(74). 

To analyze the TEP in more detail, we rewrite the expression for S by using the quasiparticle representation of the 
Green function, eq.(15), which is possible at sufficiently low temperatures in the Fermi liquid. Using the relation 
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- / ^f, (79) 

where «Sfc represents the Fermi surface and k± is the momentum perpendicular to the Fermi surface, we obtain the 
following expression: 



s _ en 2 klT 1 f dS k d f «kA | 

3ct (2tt) 3 J z k |i; k | dk±_ \ |w k |7k J e=£ , ' 



(80) 



where we performed the e-intcgration first by assuming the relation 7 <C T. In an anisotropic system, the k-dependence 
of the integrand in eq.(80) may be strong. In high-T c cuprates, for example, it is known that the anisotropy of 7k (0) 
on the Fermi surface is very large because of the strong antiferromagnctic fluctuations. The point on the Fermi surface 
where 7k takes its minimum value is called the "cold spot" , and the electrons around the cold spot mainly contribute 
to the transport phenomena. Because 7k(e£) nas a huge fcj_-dependence in high-T c cuprates around the cold spot, the 

sign of S is almost determined by the sign of e 9 \ { ^ } at the cold spot [31]. 

Next, we discuss the Nernst coefficient. Within the RTA, the Nernst coefficient is derived from eqs.(50), (29) and 
(31) by dropping all the VC's by T 22 . In an isotropic system, v by RTA is expressed in a simple form as [32,33] 



_ ir 2 klT 

vrta — — 7, 

3m 



dr(E) 
dE 



(81) 

E F 



where r(e) = l/27(e) is the energy-dependent relaxation time and Ep is the Fermi energy. According to eq.(81), v is 
determined by the energy-dependence of the relaxation time. 

Unfortunately, eq. (81) will be too simple to analyze realistic metals with (strong) anisotropy. For that purpose, 
we perform the e-integration of a xy in eq.(51) by using the quasiparticle representation. The obtained expression for 
a xy is given by 



D e 2 A 2 r 1 f ds k d _d_(.J± 

KXV 12 (2tt)3 y z k K| dk± H k 9k« [ 7k 




(82) 



where Q k ( e ) = Qk(e)/e at zero temperature. We stress that Qk( e = 0) is finite at T = as explained is §11, which 
leads to the relation v ~ 0(T7 _1 ). We stress that <2 k (0) ^ s n °t equal to Jk(0) in general, because the VC's for heat 
current and the electron one work in a different way; see discussions in §IV and §V A. 
We also comment that the Mott formula type expression for a xy , 



RTA 71 fc B J 

xy ~ 3e 



da xy (E) 



dE 



(83) 



is obtained within the RTA, by assuming that Qk(e) = eJk(e). This assumption, however, will be totally violated 
once we take the VC's into account. As a result, eq.(83) is no more valid in a correlated electron system. 

Finally, we discuss the Nernst coefficient in high-T c cuprates which increases drastically below the pseudo-gap 
temperature, T*. According to the numerical analysis based on the conserving approximation [21], fcii -dependence 

of I Jk| becomes huge due to the VC caused by the strong superconducting fluctuations. Moreover, Qk x J\t is large 
because the VC is much effective only for Jk- By considering eq.(53), the growth of the Nernst coefficient in high-T c 
cuprates under T* is caused by the enhancement of gf-|Jk|, not by J^t{E) [21]. 



VI. SUMMARY 



In the present paper, we have derived the general expressions for S, n and v in the presence of electron-electron 
interactions based on the linear response theory for the thermoelectric transport phenomena. Each expression is 
"exact" as for the most divergent term with respect to 7 _1 . The heat velocity <7k(e), which is required to calculate 
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S, k and v, is given by the Ward identity with respect to the local energy conservation law. We have studied the 
analytical properties of <7k(e) as well as the total heat current Qk(e) in detail. 

The expressions for S, k and v derived in the present paper are summarized as follows. Note that they are valid 
even if the Coulomb potential U(k) has a momentum-dependence, as discussed in Appendix C. Here, e(< 0) is the 
charge of an electron. 

(i) TEP: It is better to include the "incoherent correction", which will be important in strongly correlated systems, 
as is discussed in ref. [34]. As a result, the final expression for S is given by 

5 =^E/ *(-|{)9kx(e)[|G k (e)| 2 J kx (e)-Re{G^e)K x (e)], (84) 

where a = a xx is the electric conductivity, Wk(e) = V/j(e k + RcSk(e)), <fk(e) = e^k(e), and the total electron current 
Jk(e) is given in cq.(22). 

(ii) Thermal conductivity: In the same way, we include the incoherent correction. Then, the final expression for k 
is given by 

K = ^E /v(-^)^W[|Gk( e )| 2 Q k;c ( e )-Rc{G|( e )} (Zk;E ( e )] 

-TS 2 a, (85) 

where Qk(e) is given in eq.(30). Note that Qk(e)/e ^ Jk(e), although the Ward identity <fk(e)/e = «k(e) is rigorously 
satisfied. 

(iii) The expression for v is given by eq.(39), where a xy is given by eq.(50) or eq.(51). As for a xy (and a xy ), no 
incoherent correction exists as discussed in ref. [34] . 

These derived expressions enable us to calculate the VC's in the framework of the conserving approximation. In each 
expression, the factor 2 due to the spin degeneracy is taken into account. We note that our expression are equivalent 
to that of the relaxation time approximation (RTA), if we drop all the vertex corrections in the formulae. However, 
the RTA is dangerous because it may give unphysical results owing to the lack of conservation laws. In conclusion, the 
present work gives us the fundamental framework for the microscopic study of the thermoelectric transport phenomena 
in strongly correlated electron systems. Owing to the present work, the conserving approximation for thermoelectric 
transport coefficients becomes much practical on the basis of the Fermi liquid theory. 
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APPENDIX A: ANOTHER DERIVATION OF THE HEAT CURRENT OPERATOR, EQ.(18) 

In ref. [3], the authors derived the formula for L 12 (u>i) under the condition that electron-phonon scattering and 
the impurity scattering exist. In this appendix, for an instructive purpose, we derive eq.(19) in §11 in the case of 
the on-site Coulomb interaction by using the similar technique used in ref. [3] This fact means that the heat current 
operator in the Hubbard model can be rewritten as eq.(18). 

According to the equation of motion, the following equations are satisfied: 

[H,cl(r)} 

£ " C L + % E C i-q,<T C l+q/2,-CT C k'-q/2,-a, (Al) 
K'qa 

~ C k C ka - 2" E C k+q^c\, +q/2 _ a C k ,_ q/2 _ a . (A2) 
K'qa 



d_ 



4» 



d < ^ 
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e„ +co, 



£„ h- £„■ 

FIG. 6. The diagrammatic expression for T(ie n , ie n ' 5 ^i)- 



Using given in eq.(17) and taking eqs. (Al) and (A2) into account, we see that 
<T T j«(r)j„(0)) = E^V° (^cL(r)c feCT (r)>(0)) 

fc(7 

fcfc' gcr 

E<.( T ^LW^(r')>(0)). (A3) 

ka 

By inputting the above expression in eq.(8), we can obtain the same expression as eq.(19). As a result, j®(p = 0,w;) 
can be expressed as eq.(18). We note that eq. (18) is not exact in the case of the finite range interactions. Nonetheless, 
cq. (18) is valid for the analysis of the transport coefficient as for the most divergent term with respect to 7 -1 , as 
discussed in §111 or in Appendix C. 



APPENDIX B: DEFINITION OF T' m (pe,pV) 

Considering the convenience for readers, we list the expression for T' m (pe, pV) introduced by Eliashberg in eq.(12) 
of ref. [6], following the advice by referees. Here we dropped the momentum suffixes for simplicity. By taking the 
limit of lo — > 0, they are given by 



q~12 
q~13 

q~21 

q-23 
q~31 

<j-32 
^-33 



e, e 



e, e' 



th- 



2T 



n + cth^^-rij, 
th^ria-cth^^-ria), 



2T 

th^r 21 , 



(cth- 



^-th^^ + ccth^tf-cthi^)^ 



-th-r 23 , 



th- 



2T 



31 



cth^n-r^), 



= 0, 



-tii^r^-cth^^-r^), 



(th 



2T 



cth^tf)r 2 v 



(Bl) 



where = T^ n (e, e') (l,m — 1,2, 3, N = I, II, III, IV) is a four-point vertex function, which is introduced by the 
analytic continuation of the four-point vertex function T(ie n , ie n r, iuj{) as shown in fig. 6. For instance, r' 1 (e, e') 
comes from the analytic region [(1, 1), I] in Fig. 7 (where u>i > 0) and taking the limit ui — > at the final stage. There 
analytic properties are well studied in ref. [6]. 



APPENDIX C: THE WARD IDENTITY IN THE CASE OF THE 
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(3,1) 



\ II 



(3,2) 



II 



(3,3) 



(2,1) 



line 

4^ 



\ II / 

(2,2) ;' in 
/' IV \. 



t 



(2,3) 



(1,1) 
11/' 

/\ Ime' 



(1,2) 

-«-Im(e+i(B ; )=0 



\ I 
II \ 

(1,3) 



Im(8+iC0,)=0 

FIG. 7. The definition of the region [(l,m),N]. 



NON-LOCAL ELECTRON-ELECTRON INTERACTION 



In this appendix, we show that the expressions for S, n and v derived in the present paper is valid beyond the 
on-site Coulomb interaction. For that purpose, we reconsider the Ward identity for the following Hamiltonian H with 
a long-range interaction U(x — y) : 



h(z) = J(z) ( \^ 2 \ c(z) + ±c*(z)c(z) I vUr),{r)V(z - r)d l r. 



H 



J h{z)dz 



(CI) 
(C2) 



where V(x — y) = U(x — y)S(xo — yo), and h(z) is the local Hamiltonian. Hereafter, we drop the spin suffixes for 
simplicity. In the same reason, we put fj, = 0. 
Here, it is easy to check that 



n 2 v 2 



i 



[h(z), c(x)]5(zo — xq) = S 4 (x — z)— c(x) — 5 4 (z — x)-c(z) j c\r)c(r)V(z — r)d 4 r 

Zi 771 Z 

-±J(z)c(z)J(z)V(z-x), 



(C3) 



[H,c[x)\ 



ft 2 V 2 
2m 



c{x) — / c'(r)c(r)V(x — r)d A r ■ c(x). 



By comparing eqs.(C3) and (C4) and using the kinetic equation [H, c(x)] = —i-r£-c(x), we obtain that 

i 

[h(z),c(x)]S(zo - x ) = S 4 (z - x)(-i)-g—c(x) - -c> (z)c(z)c(x)V(z - x) 



+5 4 (z — x)i J (r)c(r)V(r — x)d 4 r ■ c(x). 



As a result, 



(T[h(z),c(x)]c\y))S(x - z ) = S 4 (x - z)-^G(x - y) - Y(x, y; z)V(x - z) 



+S 4 (x-z) Y(x,y;r)V(x-r)d\ 



(C4) 



(C5) 



(C6) 



where T is a time-ordering operator, and 



17 



Wf(k)G(k) + G(k)Wf(k) 

k+k' k k k+k' 

FIG. 8. 

The diagrammatic expression for eq.(Cll). 



/ -> k ~ A(k+k',k) A(k,k+k') , k -<, J 

w»=2k>$ A A k + , A ; \ 



FIG. 9. 

The diagrammatic expression for Aq(k) given by eq.(C13). 

Y(x,y;z) = ^(Tc\z)c(z)c(x)c\y)) 



= \JJ G(x, x')A (x', y', z)G(y', y)d i x'd 4 y', (C7) 



where A (x', y', z) is the three-point vertex function for the electron density; p(z) = c^(z)c(z). In the same way, 
(Tc(x)[h(z), c\y)])6(y - z ) = 5 4 {y - z)^-G{x - y) - Y(x, y; z)V(y - z) 

oyo 



+S A (y -z) J Y(x, y; r)V(y - r)d A r. (C8) 

As a result, we find that the following correction term 

C(x, y- z) = -Y(x, y- z)(V(x - z) + V(y - z)) 

+ j Y(x, y; r){5 A {x - z)V{x - r) + S i (y - z)V(y - r))d 4 r (C9) 

is added to eq.(34) when U(x — y) is a finite-range potential. It is easy to see that C(x,y; z) = if the potential is 
local (i.e., U(x — y) = Uq5(x — y)). 

Now, we take the Fourier transformation of C(x, y; z) according to eq.(33). If we put p^ = for [i ^ i (i=l,2 or 3) 
and \pi \ <C 1, C(k; k + pi) is given by 

C(k;k+ Pi ) = j d i xd A yd i z C(x,y;z)e lk{x ~ v)+lp{x - z) 

= %pi J d 4 xd 4 yd 4 z Y{x, y; z) [V(x - z) ■ ( Xi - z t ) + V(y - z) ■ ( Vi - Zi )] e ik ^ 

+0(p 2 t ). (CIO) 
Because the Fourier transformation of V(x) ■ x { is given by — , 

lim c(k;k+ Pi ) = r d4k , (y(fc + k) + y(fc _ k + k>)) 3 ^ 

Pi^o ipi J ok^ 

= \(w?\k)G{k) + G{k)W?\kj), (Cll) 
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which is diagrammatically shown in fig. 8. wf\k) and W^ r \k) are introduced in the last line of eq.(Cll). 

Here we note that the energy dependence of W^ ,r \k) around the Fermi level is same as that of X(fc), that is, 
ReWi(k) <~ const, and lmWi(k) « k% for \ko\ <C 1. This fact is easily recognized because 

J d A k' (Y(k + k'; k) + Y(k; k + k')) iU(k') = i (S(jfe)G(fc) + G(fc)S(fc)) , (C12) 

which is same as eq.(Cll) except for the momentum derivative on U. 

W^ r ^ given in eq.(Cll) provides the correction for the heat velocity due to the non-locality of h(z), which we 
denote as Aq(k). As shown in eq.(21) or eq.(29), in the most divergent term for L u , the (heat) current is connected 
with g( 2 \k) = |G(fc)| 2 after the analytic continuation. Bearing this fact in mind and using T 2 * = (T u + T 3l )/2 under 
the limit of T — > T as discussed in ref. [6], Aq\k) is given by 

A<f(fc) = 2G^ ({ RcW f 1 GR + G A {RcWf >}) 

= ReG _1 (fc) ■ ReWj(fc), (C13) 
which should be added to q(k) given by eq.(37) in §111. Wi(k) is given by 
W l (k) = wP(k) = w} r) {k) 

= X - J d 4 k' (A (fc + k'; k) + A (fc; k + k')) G(k + k')(-i)-^U(k'), (C14) 

which is expressed in fig. 9. In conclusion, the Ward identity for the heat velocity for general two-body interaction is 
given by <7k(e) = evk(e) + Aq(k), instead of eq.(37). 

Finally, we study the contribution of Aq(k) to the transport coefficients. Let us assume that 7 -C T at sufficiently 
lower temperatures. In this case, the correction term for the TEP due to Aq(k), AS, is calculated as 



1 W df\ z k Ag a (k,6*)J a (k, £ *) 
--o\{-Te)^ TkR) ' (C15) 

where is the quasiparticle spectrum given by the solution of RcG _1 (k, e£) = 0. Considering that Aq x (k 7 e^) = 
because of eq.(C13), we recognize that AS given by eq.(C15) is zero. 

In summary, when ?7(k) is momentum-dependent, the corrections term for the heat velocity Aq(k), given by 
eq.(C13), emerges. Fortunately, its contribution to transport coefficients would be negligible when the concept of the 
quasiparticle is meaningful, except for very high temperatures. In conclusion, the derived expressions for S, k and p, 
given by eqs.(84), (85) and (51) respectively, are valid for general electron-electron interactions, with the use of the 
heat velocity q(k) in eq.(37). 



APPENDIX D: ANOTHER PROOF OF THE WARD IDENTITY: 
BASED ON THE DIAGRAMMATIC TECHNIQUE 

In the present Appendix, we give another proof that the following generalized Ward identity is correct in a tight- 
binding model with on-site Coulomb interaction: 

e[£ k+p (e) - E k (e)] 

= T]T r J (ke; k + p, e|k' + p, e; k'e')[G k , +p (e') - G k ,(e')]e', (Dl) 

c'k' 

where T 1 is irreducible with respect to a particle-hole channel. Equation (Dl) is shown diagrammatically in fig. 10. 
Hereafter, we write that k — (e, k, a). The n-th order skeleton diagrams for the self-energy are given by 



W_ V A{P \ 

P,{ki) 

xll(k, k ai I k a2 , k a3 )U (k a4 , k a& \ k a& , fc Q7 ) • • • U {k a4n _ 4 , k a4n 3 \ ka in -2 > 

k), (D2) 



19 



e(E k+[ (e)-E k (e)) = 



k,e 

k',£' 





- k+p,£ 


k,6H 






k'+p,E' 


k',£' 
■ 




- ' 





k+p,e 

k'+p,£' 



e'Gk'(e') 

FIG. 10. The generalized Ward identity with respect to the heat velocity. 



k'-p, e; 



k+p,e k,£ 
k',£' k',£' 



e'G k .(e') 



k+p,E 

k'+p,e' 



e'Gk'(e') 



FIG. 11. 

The diagrammatic expression for eq.(D3). Here, the momentum conservation is violated by p only at the junction 
pointed by the arrow. 



where P represents the permutation of (4n — 2)-numbers, (ai, a 2 , • • • , ct4 n -2) = P (1) 1> 2, 2, • • • , In — l,2n — 1). 
A(P) = ±1 for a skclton diagram, and A{P) = for others. U(k\, k 2 \kz, k 4 ) is the two-body interaction where 
ki,k 3 are incoming and k 2l k% are outgoing, respectively [27,35]. Here we consider the on-site Coulomb interaction: 
U(k 1 ,k 2 \k 3 ,k 4 ) = US kl+ 

k3,k2+k4 ( 5ci+e3,e2+e4^o'i,tT2^(T3,<T4^ai,— era- Note that in eq.(D2), the tadpole-type (Hartree-type) 
diagrams are dropped because they are fc-independent. 

Next, we consider the right-hand-site of eq.(Dl), which can be rewritten as 



tJ2 [r 7 (fc; k + P \k', k'-p)- r 7 (fc ; k + p\k' + P , k')] G v {<?y, 



(D3) 



e'k' 



which is shown in fig.ll. 

Then, the n-th order skeleton diagrams for eq.(D3) is expressed as 

^ ~j '*G , fe 1 Gfe 2 • • • Gk 2n _ 1 



P,{ki} 
X 



^ ^aif^(fc,^ai p|^«2 > k a3 ) + c ai U(k, k ai \k a2 k a3 p) e a2 U{k, k ai \k a2 + p, k a3 ) 
xU(k a4 , k a5 \k a6 , k aj ) ■ ■ ■ U(k a4n4 , k a4n3 \k a4n _ 2 , k + p) 

£c«5 U (ka 4 , p| ^q 6 , ^a 7 ) ~\~ €aj U {ka. 4 , /Ca5 | ka 6 , fc ft7 p) 

~^aiU (k a4 + p, fc as |fc a6 , k ar ) — e ae U(k a4 , k a5 \ k ae + p, k ar ) 
xll(k, k ai \k a2 , k a3 ) • • • U(k a4n _ 4 , k a4n3 \k a4n _ 2 , fc + p) 



+ ■ 

+ 



^a4 n _3 U{k a4n 4 , k a4n 3 P|^a4„-2 7 k "T" p) e ft4 )l _4^(^a4 Tl _4 ~t" P, k a4n _ 3 | k a4n _ 2 , 
— e a 4 „_2^(^a 4 „_4 , ^Ct4n-3 l^"4n-2 + Pi ^ + P) 

xll(k, k ai \ k a2 , /c Q3 ) • • • U{k a4n _ s , k a4n 7 \ k a4n _ 6 , fc a4n _ 5 ) > . 



fc + p) 



(D4) 



Because of the relation £/ (fci + p, A^lfe, fct) = U(k\, k 2 \k 3 + p, k 4 ) = U(ki,k 2 — p\k%, k 4 ) — U(k\, k 2 \k 3 , k 4 — p) and 
the energy conservation with respect to U, we see that eq.(D4) + e(Sfe — £fc+p) = 0. As a result, the generalized 
Ward identity, eq.(Dl) is proved in the framework of the microscopic perturbation theory. 

By taking the limit |p| — > of eq.(Dl), we obtain 
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e„ • V fc S k (e„) = T ^ lim r J (ke„; k + p, e n \k' + p, e n r, k'e„/) • Gk'+p(e«') G k' 

xe n / (t£, + Vfc/£ k ' . (D5) 
Equation (D5) is rewritten by using the reducible four-point vertex T as 

e„ • V fc Sk(e„) = TV lim T(ke„; k + p, e„|k' + p, e n /;k'e„') • Gk'+p(en') G k'(e«') • e«'^- (D6) 

e n /k' 

After the analytic continuation of eq.(D6), we find that 

k',i=l,3 



where the four-point vertex T y is introduced by Eliashberg in ref. [6], and ^(e) is introduced in §11. Thus, the Ward 
identity cq.(D7), is derived in terms of the diagrammatic technique. Equation (D7) is equivalent to eq.(37), which is 
expressed in terms of the zero temperature perturbation method. 

We note that "the usual Ward identity" related to the charge conservation law is given by [26,27] 

V fe S^( e )= J2 /^T»(kc,kV)(/«(cX'- (D8) 

k', i=l,3 

According to eqs. (D7) and (D8), 

£ /^T«(ke,kV) ff «( e ')-^e' = e- ^ /^(ke,kV)^V)^- (D9) 

k',i= 1,3 k,i— 1,3 
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